Дискуссионный математический форумМатематический форум

Математический форум Math Help Planet

Обсуждение и решение задач по математике, физике, химии, экономике

Теоретический раздел
Часовой пояс: UTC + 4 часа [ Летнее время ]
MathHelpPlanet.com RSS-лента Математического форума

Часовой пояс: UTC + 4 часа [ Летнее время ]


Z-преобразование и его свойства

Z-преобразование и его свойства


Основные определения


1. Оригинал — последовательность

{f(k), k=0,1,…}
, удовлетворяющая условию:
|f(k)|<Meσk
, где
M
и
σ
— положительные постоянные (рис. 5.22).


2. Изображение последовательности

{f(k), k=0,1,…}
— функция
F(z)
комплексного переменного
z
, определяемая равенством


F(z)=∑k=0∞f(k)zk.(5.58)

Изображение является аналитической функцией при

|z|>eσ
.


Совокупность всех оригиналов называется пространством оригиналов, а совокупность всех изображений — пространством изображений.


3. Переход, определяющий изображение

F(z)
по оригиналу
{f(k), k=0,1,…}
, называется Z-преобразованием:


F(z)=Z[f(k)].
(5.59)

4. Оригинал по изображению находится с помощью обратного Z-преобразования по формуле:


f(k)=12πi∮CF(z)zk−1dz,k=0,1,…,
(5.60)

где

C
— контур, внутри которого лежат все особые точки функции
F(z)
.


Замечания 5.14


1. Название Z-преобразование определяется буквой

z
, выбранной для обозначения переменной. Такое название противоречит существующему обычаю называть часто применяемые преобразования по имени ученого. В некоторых источниках Z-преобразование называется преобразованием Лорана, так как ряд (5.58) дает разложение функции
F(z)
в ряд Лорана в окрестности бесконечно удаленной точки.


2. Z-преобразование можно рассматривать как частный случай преобразования Лапласа, а именно как преобразование в пространстве ступенчатых оригиналов. Ступенчатая функция — это разрывная функция целочисленного аргумента, которая в общем случае имеет разрывы при каждом натуральном значении аргумента, оставаясь между ними постоянной. По рис. 5.22 можно представить такую функцию

f(t)
, для которой
f(t)=f(0)
при
t∈[0;1), f(t)=f(1)
при
t∈[1;2)
и т.д. Оригиналом по Лапласу такая функция является, если выполняется условие
|f(t)|<Meσt
, то есть
|f(k)|<Meσk
. Ее можно представить в виде ряда
f(t)=∑k=0∞f(k)[1(t−k)−1(t−(k+1))]
. Тогда по теореме запаздывания


L[f(t)]=∑k=0∞f(k)L{[1(t−k)−1(t−(k+1))]}=∑k=0∞f(k)[e−kpp−e−(k+1)pp]=1−e−pp∑k=0∞f(k)e−kp,

т.е. нахождение преобразования Лапласа в рассматриваемом классе оригиналов сводится к нахождению суммы ряда Тейлора:


Φ(u)=∑k=0∞akuk,ak=f(k),u=e−p.

Обозначая

z=ep
(тогда
z=1u
) и отбрасывая множитель
1−e−pp
, общий для всех оригиналов данного класса, получаем соотношение (5.58):
F(z)=Φ(1z)=∑k=0∞f(k)zk
. При этом обратное преобразование — есть задача нахождения коэффициентов' разложения функции
Φ(u)
в ряд Тейлора по степеням
u
. Коэффициенты вычисляются по формуле (3.16):
ak=12πi∮CΦ(u)uk+1du, k=0,1,…
или (3.17):
ak=Φ(k)(0)k!
. При
u=1z
получаем ряд Лорана (5.58) функции
F(z)=Φ(u)
в окрестности
z=∞
с коэффициентами (5.60):
f(k)=12πi∮CF(z)zk−1dz, k=0,1,…
(здесь учтено, что
du=−1z2dz
, направление обхода контура при преобразовании
z=1u
меняется на противоположное и интеграл меняет знак при изменении направления обхода).




Свойства Z-преобразования


Положим, что

F(z)=Z[f(k)]
.


1. Линейность. Для любых постоянных

ci, i=1,2,…,m
справедливо


Z[c1f1(k)+…+cmfm(k)]=c1F1(k)+…+cmFm(k),
(5.61)

где

{f1(k),k=0,1,…},…,{fm(k),k=0,1,…}
— оригиналы, a
F1(z),…,Fm(z)
— их изображения.


2. Запаздывание (формула запаздывания), где

f(k−n)=0
при
k−n<0
(рис. 5.23):


Z[f(k−1)]=z−1F(z),
(5.62)

Z[f(k−n)]=z−nF(z), n=1,2,…,
(5.63)

3. Опережение (формула опережения):


Z[f(k+1)]=z[F(z)−f(0)],
(5.64)

Z[f(k+n)]=znF(z)−znf(0)−…−zf(n−1).
(5.65)

4. Дифференцирование изображения:


Z[k⋅f(k)]=−z⋅dFdz.
(5.66)

5. Умножение изображений. Свертке оригиналов соответствует произведение изображений:


Z[∑i=0kf(i)⋅g(k−i)]=F(z)⋅G(z),
(5.67)

где

G(z)=Z[g(k)]
. Сверткой оригиналов
{f(k),k=0,1,…}
и
{g(k),k=0,1,…}
называется сумма


∑i=0kf(i)⋅g(k−i).
(5.68)

6. Теоремы о предельных значениях. Если

F(z)=Z[f(n)]
, то


f(0)=limz→∞F(z),f(1)=limz→∞z[F(z)−f(0)],f(2)=limz→∞z2[F(z)−f(0)−f(1)⋅z−1],⋮
(5.69)

где

z
стремится к бесконечности вдоль произвольного пути. Если
limk→∞f(k)
существует, то


limk→∞f(k)=limz→1(z−1)F(z).
(5.70)

Заметим, что первую формулу в (5.69) можно получить из (5.58):


F(z)=f(0)+f(1)z+f(2)z2+…+f(k)zk+…
при
z→∞
.

Умножая обе части

F(z)
на
z
, получаем
f(1)=limz→∞z[F(z)−f(0)]
и т.д.




Нахождение изображения по оригиналу


Для нахождения изображений используется (5.58) и свойству Z-преобразования.


Пример 5.61. Найти изображения по оригиналам: а)

f(k)=ak
; б)
f(k)=1
; в)
f(k)=akk!
; г)
f(k)=eαk
.


▼ Решение

а) По формуле (5.58) имеем


F(z)=∑k=0∞akzk=11−a⧸|z=zz−a,

где использована формула суммы бесконечно убывающей геометрической прогрессии с первым членом

a1=1
и знаменателем
q=az
при
|q|<1
, то есть
|z|>|a|
.


б) Аналогично пункту "а"


F(z)=∑k=0∞1zk=11−1⧸|z=zz−1.

в)

F(z)=∑k=0∞akk!zk=ea⧸|z
при
|z|>|a|
; где использовано типовое разложение.


г)

F(z)=∑k=0∞eαkzk=11−eα⧸|z=zz−eα
при
|z|>|eα|=eRe⁡α
. Заметим, что при
α=0
имеем результат, полученный в пункте "б".


В табл. 5.2 приведены наиболее часто встречающиеся в примерах соответствия при Z-преобразовании. Формулы 1, 7, 11 получены в примере 5.61.


Таблица 5.2. Таблица Z-преобразований


№
f(k)
F(z)
№
f(k)
F(z)
1
1
zz−1
11
akk!
ea⧸|z
2
(−1)k
zz+1
12
k(k−1)2
z(z−1)3
3
k
z(z−1)2
13
1m!∏i=0m−1(k−i)
z(z−1)m+1
4
k2
z(z+1)(z−1)3
14
sin⁡kβ
zsin⁡βz2−2zcos⁡β+1
5
kak−1
z(z−a)2
15
cos⁡kβ
z(z−cos⁡β)z2−2zcos⁡β+1
6
Ckmak
amz(z−a)m+1
16
aksin⁡kβ
azsin⁡βz2−2azcos⁡β+a2
7
ak
zz−a
17
akcos⁡kβ
z(z−acos⁡β)z2−2azcos⁡β+a2
8
(k+1)ak
z2(z−a)2
18
−aksin⁡(k−1)β
zsin⁡β(z−2acos⁡β)z2−2azcos⁡β+a2
9
{0,k=01⧸|k,k≠0
ln⁡1z−1
19
aksh⁡kβ
azsh⁡βz2−2azch⁡β+a2
10
{0,k=0(−1)k⧸|k,k≠0
ln⁡z+1z
20
akch⁡kβ
z(z−ach⁡β)z2−2azch⁡β+a2



Нахождение оригинала по изображению


Для нахождения оригинала используются следующие способы:


1) разложение функции
F(z)
в ряд Лорана в окрестности бесконечно удаленной точки;
2) формула (5.60);
3) табл. 5.2 и свойства Z-преобразования;
4) формула (5.69);
5) формула
f(k)=1k!(dkF(z−1)dzk)|z=0,k=0,1,2,…
(5.71)

6) если

F(z)
представляет собой рациональную функцию, т.е. отношение двух многочленов:
F(z)=P(z)Q(z)
, где степень многочлена в силу аналитичности
F(z)
при
z=∞
не превышает степени многочлена
Q(z)
, то можно разделить многочлен
P(z)
на
Q(z)
одним из обычных способов. Это хотя и не даст общего выражения для
f(k)
, но позволит численным путем определить сколь угодно большое число значений
f(k)
.


Замечания 5.15


1. Сравнивая (5.60) с определением (4.16) вычета в

z=∞
, получаем (где
zm
— особые точки функции
F(z)
)


f(k)=−resz=∞⁡F(z)zk−1,f(k)=∑m=1nresz=zm⁡F(z)zk−1.

2. Формула (5.71) получается из формулы для коэффициентов ряда Тейлора

?k=1k!Φ(k)(0)
(см. замечание 5.14), где
Φ(1z)=F(z), Φ(0)=F(1z)z=0
.


Пример 5.62. Найти оригиналы для функций:
а)
F(z)=z(z−1)2
; б)
F(z)=z2z2−z−1
; в)
F(z)=1(z−a)(z−b)
;
г)
F(z)=z(z2−9)2
; д)
F(z)=z2+zz2−z+1
; е)
F(z)=z2z2−2z+1
;

▼ Решение

а) Первый способ. Разложим функцию

F(z)
в ряд Лорана в окрестности бесконечно удаленной точки, учитывая, что
1(z−1)2=(11−z)′:


F(z)=zddz(11−z)=zddz(−1z(1−1⧸|z))=−zddz[1z(1+1z+1z2+…)]==−zddz(1z+1z2+1z3+…)=−z(−1z2−2z3−3z4−…)==1z+2z2+3z3+…=∑k=1∞kzk.

Сравнивая с (5.58), имеем

f(k)=k, k=1,2,…
По формуле (5.69) получаем
f(0)=limz→∞F(z)=0
. Поэтому
f(k)=k, k=0,1,2,…
.


Второй способ. По формуле (5.60) получаем


f(k)=12πi∮Cz(z−1)2zk−1dz,k=0,1,2,…

В точке

z=1
— полюс второго порядка; интеграл можно брать по любому контуру
C:|z|=R
, где
R>1
, так как в
|z|<R, R>1
содержится единственная особая точка функции
F(z)
. По формуле (4.19)


∮Czk(z−1)2dz=2πiresz=1⁡zk(z−1)2,

где согласно (4.22)

resz=1⁡zk(z−1)2=limz→1ddzzk=limz→1kzk−1=k
. Поэтому


f(k)=12πi⋅2πik=k
, то есть
f(k)=k, k=0,1,2,…
.

Заметим, что можно было использовать п.1 замечаний 5.15:


f(k)=resz=1⁡F(z)zk−1=resz=1⁡zk(z−1)2=lim(zk)′=k.

Третий способ. Представим

F(z)
в виде рациональной функции:
F(z)=zz2−2z+1=P(z)Q(z)
. Поделим числитель на знаменатель, в результате получаем


F(z)=1z+2z2+3z3+…
и
f(k)=k, k=0,1,2,…

Четвертый способ. По формуле (5.69) имеем


f(0)=limz→∞zz2−2z+1=0,f(1)=limz→∞z(zz2−2z+1−0)=1,f(2)=limz→∞z2(zz2−2z+1−0−1⋅1z)=limz→∞z4−z2(z2−2z+1)(z2−2z+1)z=limz→∞2z3−z2z3−2z2+z=2,

и так далее:

f(k)=k, k=0,1,2,…
.


б) По формуле (5.60) имеем


f(k)=12πi∮Cz2z2−z−1zk−1dz,0,1,2,…

Особые точки

z1=1+52, z2=1−52
полюсы первого порядка (корни уравнения
z2−z−1=0
). По формуле (4.19) имеем


f(k)=12πi⋅2πi[resz=z1⁡zk+1(z−z1)(z−z2)+resz=z2⁡zk+1(z−z1)(z−z2)].

По формуле (4.23) для вычисления вычета в полюсе первого порядка получаем


f(k)=resz=z1⁡zk+1(z−z1)(z−z2)+resz=z2⁡zk+1(z−z1)(z−z2)==limz→z1(z−z1)zk+1(z−z1)(z−z2)+limz→z2(z−z2)zk+1(z−z1)(z−z2)==z1k+1−z2k+1z1−z2=15[(1+52)k+1−(1−52)k+1].

в) Представим

F(z)
в виде суммы двух элементарных дробей и разложим каждую из них в окрестности
z=∞:


F(z)=1(z−a)(z−b)=1a−b(1z−a−1z−b)=1(a−b)z(11−a⧸|z−11−b⧸|z)==1(a−b)z∑k=0∞(akzk−bkzk)=1a−b∑k=0∞ak−bkzk+1=1a−b∑k=2∞ak−1−bk−1zk.

Поэтому

f(k)=ak−1−bk−1a−b
при
k⩾2
. По формуле (5.69) находим
f(0)=limz→∞F(z)=0,
f(1)=limz→∞zF(z)=0
. Окончательный ответ:
f(k)=ak−1−bk−1a−b
при
k⩾1
.


г) Разложим функцию

F(z)⧸|z
на простейшие дроби:


F(z)z=1(z2−9)2=136[1(z−3)2+1(z+3)2]−1108(1z−3−1z+3),

поэтому

F(z)=1108[3z(z−3)2+3z(z+3)2−zz−3+zz+3].

По формулам 5,7 из табл. 5.2 и


Z−1[3z(z−3)2]=3k3k−1=k3k;Z−1[3z(z+3)2]=3k(−3)k−1=−k(−3)k;Z−1[zz−3]=3k;Z−1[zz+3]=(−3)k.

Используя свойство линейности, имеем


f(k)=1108[k3k−k(−3)k−3k+(−3)k]=1108(k−1)3k[1k−(−1)k],  k=0,1,2,…

д) Заметим, что дискриминант трехчлена в знаменателе отрицателен. Наиболее близкими по форме записи к данному изображению являются функции, соответствующие формулам 14 и 15 из табл. 5.2. Чтобы выделить значения

sin⁡β
и
cos⁡β
, представим изображение в виде


F(z)=z2+zz2−z+1=z(z−12)+32zz2−2z12+1=z(z−12)z2−2z12+1+32zz2−2z12+1,

При

cos⁡β=12
, то есть
β=π3
, первое слагаемое совпадает с изображением в формуле 15 из табл. 5.2. Так как
sin⁡β=sin⁡π3=22
, то преобразуем второе слагаемое к табличному. Получим


F(z)=z(z−cos⁡π3)z2−2zcos⁡π3+1+3⋅zsin⁡π3z2−2zcos⁡π3+1.

По формулам 15, 14 из табл. 5.2 и по свойству линейности


f(k)=cos⁡πk3+3sin⁡πk3,k=0,1,2,…

е) Представим изображение в виде


F(z)=z2z2−2z+1=z⋅zz2−2z22+1=2z⋅zsin⁡π4z2−2zcos⁡π4+1.

По формуле 16 из табл.5.2 и свойству линейности


Z−1[F(z)z]=Z−1[2⋅zsin⁡π4z2−2zcos⁡π4+1]=2sin⁡πk4=φ(k).

По формуле опережения (5.64) с учетом

φ(0)=0
получаем:


f(k)=Z−1[F(z)]=2sin⁡(k+1)π4.



Области применения Z-преобразования


Решение линейных разностных уравнений с постоянными коэффициентами
Пусть заданы:
а) линейное разностное уравнение:

anx(k+n)+an−1x(k+n−1)+…+a1x(k+1)+a0x(k)=g(k),k=0,1,2,…
(5.72)

где

anan−1,…,a1,a0
— постоянные коэффициенты;
n
— целое число; последовательность
{g(k),k=0,1,2,…}
;


б) начальные условия:


x(0)=x0, x(1)=x1, …, x(n−1)=xn−1.
(5.73)

Требуется найти решение

x(k)
разностного уравнения, удовлетворяющее начальным условиям.


3амечание 5.16. Аналогично формируется постановка задачи решения систем линейных разностных уравнений.


Будем предполагать, что последовательности

{g(k),k=0,1,2,…}
и
{x(k),k=0,1,2,…}
принадлежат пространству оригиналов. Для решения поставленной задачи можно применить аппарат операционного исчисления — метода решения задач, содержащего следующие этапы, отраженные на рис. 5.24.




Алгоритм решения задачи (5.72),(5.73)


1. От известных и неизвестных последовательностей (оригиналов) перейти к их Z-изображениям. Записать уравнение (систему) в изображениях, соответствующее решаемой задаче.

2. Решить полученное уравнение (систему): найти изображение искомого решения.

3. Применить обратное Z-преобразование: найти оригинал для полученного в п.2 изображения.


▼ Примеры 5.63-5.71

Пример 5.63. Найти решение задачи:

x(k+2)−x(k+1)+x(k)=0, x(0)=1, x(1)=2
.


▼ Решение

1. Перейдем от оригиналов к изображениям, используя формулы (5.61), (5.64), (5.65):


Z[x(k)]=X(z);Z[x(k+1)]=zX(z)−zx(0)=zX(z)−z;Z[x(k+2)]=z2Z(z)−z2x(0)−zx(1)=z2X(z)−z2−2z;Z[0]=0.

Запишем уравнение для изображений:

z2X(z)−z2−2z−zX(z)+z+X(z)=0
.


2. Решим уравнение для изображений:


(z2−z+1)X(z)=z2+z⇒X(z)=z2+zz2−z+1.

3. Найдем оригинал для

X(z):


X(z)=z(z−12)+32zz2−2z12+1=z(z−cos⁡π3)z2−2zcos⁡π3+1+3⋅zsin⁡π3z2−2zcos⁡π3+1.

По формулам 14 и 15 из табл. 5.2 получаем (см. пример 5.62, пункт "д")


x(k)=cos⁡πk3+3sin⁡πk3,k=0,1,2,…

Легко проверить, что при

k=0
и
k=1
начальные условия выполняются.


Пример 5.64. Найти решение задачи:

x(k+2)−3x(k+1)+x(k)=0, x(0)=12, x(1)=32
.


▼ Решение

1. Перейдем от оригиналов к изображениям, используя формулы (5.61), (5.64), (5.65):


Z[x(k)]=X(z);Z[x(k+1)]=zX(z)−zx(0)=zX(z)−z2;Z[x(k+2)]=z2X(z)−z2x(0)−zx(1)=z2X(z)−z22−z32;Z[0]=0.

Запишем уравнение для изображений:


z2X(z)−z22−z32−3zX(z)+32z+X(z)=0.

2. Решим уравнение для изображений:


(z2−3z+1)X(z)=z22⇒X(z)=z22(z2−3z+1).

3. Найдем оригинал для

X(z):


X(z)=z22(z2−2z32+1)=z⋅zsin⁡π6z2−2zcos⁡π6+1.

По формуле 14 из табл. 5.2 и (5.64) получаем (см. пример 5.62, пункт "е")


x(k)=sin⁡π(k+1)6,k=0,1,2,…

При использовании формулы (5.64) было учтено, что


Z−1[X(z)z]=Z−1[zsin⁡π6z2−2zcos⁡π6+1]=sin⁡πk6=f(k),f(0)=0.

Пример 5.65. Найти решение задачи:

x(k+2)−3x(k+1)−10x(k)=0, x(0)=3, x(1)=−1
.


▼ Решение

1. Перейдем от оригиналов к изображениям, используя формулы (5.61), (5.64), (5.65):


Z[x(k)]=X(z);Z[x(k+1)]=zX(z)−zx(0)=zX(z)−3z;Z[x(k+2)]=z2X(z)−z2x(0)−zx(1)=z2X(z)−3z2+z;Z[0]=0.

Запишем уравнение для изображений:


z2X(z)−3z2+z−3zX(z)+9z−10X(z)=0.

2. Решим уравнение для изображений:


(z2−3z−10)X(z)=3z2−10z⇒X(z)=3z2−10zz2−3z−10.

3. Найдем оригинал для

X(z): X(z)=3z2−10z(z−5)(z+2)
.


Первый способ. Запишем

X(z)
в виде произведения, удобного для применения формулы 7 из табл. 5.2, и разложим дробь на элементарные:


X(z)=3z2−10z(z−5)(z+2)=z⋅3z−10(z−5)(z+2)=z()=57⋅zz−5+167⋅zz+2.

По формуле 7 из табл. 5.2 получаем


x(k)=57⋅5k+1617⋅(−2)k=175k+1+167(−2)k.

Второй способ. Согласно п.1 замечаний 5.15


x(k)=resz=5⁡3z2−10z(z−5)(z+2)zk−1+resz=−2⁡3z2−10z(z−5)(z+2)zk−1==(3z2−10z)(z−5)(z−5)(z+2)zk−1|z=5+(3z2−10z)(z+2)(z−5)(z+2)zk−1|z=−2==75−5075k−1−327(−2)k−1=175k+1+167(−2)k.

Пример 5.66. Найти решение задачи:

x(k+2)+4x(k+1)+3x(k)=1, x(0)=1, x(1)=1
.


▼ Решение

1. Перейдем от оригиналов к изображениям, используя формулы 1 из табл. 5.2, (5.61), (5.64), (5.65):


Z[x(k)]=X(z);Z[x(k+1)]=zX(z)−zx(0)=zX(z)−z;Z[x(k+2)]=z2X(z)−z2x(0)−zx(1)=z2X(z)−z2−z;Z[1]=zz−1.

Запишем уравнение для изображений:


z2X(z)−z2−z+4zX(z)−4z+3X(z)=zz−1.

2. Решим уравнение для изображений:


X(z)=z2+5z+zz−1z2+4z+3=z3+4z2−4z(z−1)(z+1)(z+3).

3. Найдем оригинал для

X(z)
. Как и в предыдущем примере, представим
X(z)z
в виде суммы простейших дробей:


X(z)z=z2+4z−4(z−1)(z+1)(z+3)=18⋅1z−1+74⋅1z+1−78⋅1z+3.

Отсюда

X(z)=18⋅zz−1+74⋅zz+1−78⋅zz+3.

По формулам 1,2 и 7 из табл. 5.2 имеем:
x(k)=18+74(−1)k−78(−3)k
.

Пример 5.67. Найти решение задачи:

x(k+2)−5x(k+1)+6x(k)=2⋅4k, x(0)=2, x(1)=1
.


▼ Решение

1. Перейдем от оригиналов к изображениям, используя формулы 7 из табл. 5.2, (5.61), (5.64), (5.65):


Z[x(k)]=X(z);Z[x(k+1)]=zX(z)−zx(0)=zX(z)−2z;Z[x(k+2)]=z2X(z)−z2x(0)−zx(1)=z2X(z)−2z2−z;Z[2⋅4k]=2⋅zz−4.

Запишем уравнение для изображений:


z2X(z)−2z2−z−5zX(z)+10z+6X(z)=2⋅zz−4.

2. Решим уравнение для изображений:


X(z)=2⋅z(z−4)(z2−5z+6)+2z2−9zz2−5z+6.

3. Найдем оригинал для

X(z)=z⋅(2(z−4)(z−3)(z−2)+2z−9(z−3)(z−2))
. Разложим каждую из дробей на элементарные:


2(z−4)(z−3)(z−2)=1z−4+1z−2−2z−32z−9(z−3)(z−2)=−3z−3+5z−2.

Тогда

X(z)z=1z−4+6z−2−5z−3⇒X(z)=zz−4+6zz−2−5zz−3.

По формуле 7 из табл. 5.2 имеем:

x(k)=4k+6⋅2k−5⋅3k
.


Пример 5.68. Найти решение задачи:

x(k+2)−4x(k+1)+4x(k)=3k, x(0)=x0, x(1)=x1
.


▼ Решение

1. Перейдем от оригиналов к изображениям, используя формулы 7 из табл. 5.2, (5.61), (5.64), (5.65):


Z[x(k)]=X(z);Z[x(k+1)]=zX(z)−zx0;Z[x(k+2)]=z2X(z)−z2x0−zx1;Z[3k]=zz−3.

Запишем уравнение для изображений:


z2X(z)−z2x0−zx1−4zX(z)+4zx0+4X(z)=zz−3.

2. Решим уравнение для изображений:


X(z)=z(z−3)(z2−4z+4)+z2x0+zx1−4x0zz2−4z+4.

3. Найдем оригинал для

X(z)
. Используем разложение


1(z−3)(z2−4z+4)=1z−3−1z−2−1(z−2)2.

Тогда

X(z)=zz−3−zz−2−z(z−2)2+x0⋅z2(z−2)2+(x1−4x0)z(z−2)2.

По формулам 5, 7, 8 из табл. 5.2 получаем


x(k)=3k−2k−k⋅2k−1+x0(k+1)2k+(x1−4x0)k2k−1==3k−2k+x0(k+1)2k+(x1−4x0−1)k⋅2k−1.

Пример 5.69. Найти решение задачи:

x(k+3)−3x(k+2)+3x(k+1)−x(k)=2k, x(0)=0, x(1)=0, x(2)=1
.


▼ Решение

1. Перейдем от оригиналов к изображениям, используя формулы 7 из табл. 5.2, (5.61), (5.64), (5.65):


Z[x(k)]=X(z);Z[x(k+1)]=zX(z)−zx(0)=zX(z);Z[x(k+2)]=z2X(z)−z2x(0)−zx(1)=z2X(z);Z[2k]=zz−2.Z[x(k+3)]=z3X(z)−z3x(0)−z2x(1)−zx(2)=z3X(z)−z.

Запишем уравнение для изображений:


z3X(z)−z−3z2X(z)+3zX(z)−X(z)=zz−2.

2. Решим уравнение для изображений:


X(z)=z+zz−2z3−3z2+3z−1=z(z−1)3+z(z−2)(z−1)3.

3. Найдем оригинал для

X(z)
. Используем представление


1(z−2)(z−1)3=1z−2−1(z−1)3−1z−1−1(z−1)2.

Тогда

X(z)=z(z−1)3+zz−2−z(z−1)3−zz−1−z(z−1)2=zz−2−zz−1−z(z−1)2.

По формулам 1, 3 и 7 из табл. 5.2 получаем

x(k)=2k−1−k
.


Пример 5.70. Найти решение задачи:

{x(k+2)−y(k)=0,y(k+2)+x(k)=0;x(0)=y(0)=1, x(1)=2, y(1)=0
.


▼ Решение

1. Перейдем от оригиналов к изображениям, используя формулы (5.64), (5.65):


Z[x(k)]=X(z);Z[y(k)]=Y(z);Z[x(k+1)]=zX(z)−zx(0)=zX(z)−z;Z[y(k+1)]=zY(z)−zy(0)=zY(z)−z;Z[x(k+2)]=z2X(z)−z2x(0)−zx(1)=z2X(z)−z2−2z;Z[y(k+2)]=z2Y(z)−z2y(0)−zy(1)=z2Y(z)−z2;z[0]=0.

Запишем систему уравнений для изображений:


{z2X(z)−z2−2z−Y(z)=0,z2Y(z)−z2+X(z)=0.

2. Решим систему уравнений для изображений. Выражая

Y(z)
из первого уравнения и подставляя во второе, получаем


z4X(z)−z4−2z3−z2+X(z)=0.

Отсюда, учитывая, что

z4+1=(z2+1)2−2z2=(z2−2z+1)(z2+2z+1)
, имеем


X(z)=z4+2z3+z2z4+1=z2(z2+2z+1)z4+1=z2z2−2z+1.Y(z)=z4z2−2z+1−z2−2z=z2−2zz2−2z+1.

3. Найдем оригиналы для

X(z)
и
Y(z)
. Как в примере 5.62, пункт "д" дискриминант знаменателя отрицателен, поэтому преобразуем дроби к виду, удобному для применения формул 14-17 из табл. 5.2:


X(z)=z2z2−2z+1=2z⋅zsin⁡π4z2−2zcos⁡π4+1,Y(z)=z2−2zz2−2z+1=z(z−22)−z22z2−2zcos⁡π4+1=z(z−cos⁡π4)z2−2zcos⁡π4+1−sin⁡π4z2−2zcos⁡π4+1.

По формулам 14, 15 из табл. 5.2 и (5.64) получаем


x(k)=2sin⁡π(k+1)4,y(k)=cos⁡πk4−sin⁡πk4.

Здесь было учтено, что


Z−1[zsin⁡π4z2−2zcos⁡π4+1]=sin⁡πk4=f(k),k=1,2,…;f(0)=0.

Пример 5.71. Найти решение задачи:

{x(k+1)−x(k)+y(k)=3k,y(k+1)+2x(k)=−3k;x(0)=3, y(0)=0
.


▼ Решение

1. Перейдем от оригиналов к изображениям, используя формулу 7 из табл. 5.2, (5.64):


Z[x(k)]=X(z);Z[y(k)]=Y(z);Z[x(k+1)]=zX(z)−zx(0)=zX(z)−3z;Z[y(k+1)]=zY(z)−zy(0)=zY(z);Z[3k]=zz−3.

Запишем систему уравнений для изображений:


{zX(z)−3z−X(z)+Y(z)=zz−3;zY(z)+2X(z)=−zz−3.

2. Решим систему уравнений для изображений. Умножая первое уравнение на

z
и вычитая второе, получаем


z2X(z)−3z2−zX(z)−2X(z)=z2z−3+zz−3.

Отсюда

{X(z)=3z3−8z2+z(z−3)(z2−z−2)=3z3−8z2+z(z−3)(z−2)(z+1)Y(z)=−1z−3−6z2−16z+2(z−3)(z−2)(z+1).

3. Найдем оригиналы для

X(z),Y(z)
. Используем представления


X(z)z=3z2−8z+1(z−3)(z−2)(z+1)=1z−3+1z−2+1z+1,6z2−16z+2(z−3)(z−2)(z+1)=2z−2+2z+1+2z−3.

Тогда

X(z)=zz−3+zz−2+zz+1,Y(z)=−2z−2−2z+1−2z−3=z−1(−2zz−2−2zz+1−2zz−3).

По формулам (5.62) и 7 из табл. 5.2 получаем


x(k)=3k+2k+(−1)k,y(k)=(−2)⋅2k−1−2⋅(−1)k−1−3⋅3k−1, k⩾1,y(0)=0⇕y(k)=−2k+2(−1)k−3k,k=0,1,2,…

Часовой пояс: UTC + 4 часа [ Летнее время ]


Яндекс.Метрика

Copyright © 2010-2016 MathHelpPlanet.com. All rights reserved